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Abstract 

We investigate the problem of extracting a static potential between a quark and its antiquark in a quark-gluon plasma 
(QGP) from lattice-QCD computations of the singlet free energy, FQQ(r). We utilize the thermodynamic T -matrix 
formalism to calculate the free energy from an underlying potential ansatz resummed in ladder approximation. Imag¬ 
inary parts of both QQ potential-type and single-quark selfenergies are included as estimated from earlier results of 
the T -matrix approach. We find that the imaginary parts, and in particular their (low-) energy dependence, induce 
marked deviations of the (real part of the) potential from the calculated free energy. When fitting lattice results of the 
latter, the extracted potential is characterized by significant long-range contributions from remnants of the confining 
force. We briefly discuss consequences of this feature for the heavy-quark transport coefficient in the QGR 

Keywords: Heavy-quark potential, Bethe-Salpeter Equation, Quark-Gluon Plasma 
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1. Introduction 


Heavy-flavor (HE) particles, containing charm or 
bottom quarks, play a central role in the analysis of 
strongly interacting matter as produced in ultrakela- 
tivistic heavy-ion collisions (URHICs) 111]. The nu¬ 
clear modification factor and elliptic flow of D and B 
mesons (as well as their semileptonic decay leptons) 
are key measures to establish and quantiW the strong¬ 
coupling nature of the QCD medium yj-ldt], while the 
production patterns of charmonia and bottomonia 
are believed to encode information on its deconflnement 
properties. Eor both observables, a good understanding 
of the in-medium interactions of the heavy quarks and 
quarkonia is mandatory. 

Lattice-QCD (IQCD) computations allow to study the 
properties of heavy-quark (HQ) systems in the quark- 
gluon plasma (QGP) from first principle. In particular, 
euclidean correlation functions of quarkonia, HQ sus¬ 
ceptibilities and static QQ free energies have been com¬ 
puted with good precision SESi. A non-perturbative 
T-matrix approach (“Brueckner theory”) for quarko¬ 
nia and heavy quarks in the QGP has been developed 
to interpret the IQCD results jlT, 12], and the calcu¬ 
lated spectral functions and transport coefficients have 
been applied to both hidden ll3n and open HP observ¬ 
ables il4ll in URHICs. The main theoretical uncertainty 
in this and other approaches based on the Schrodinger 


equation Bfil] remains the definition of the driving 
kernel or potential. Thus far, it has been bracketed by 
using either the temperature dependent free or inter¬ 
nal energies (or combinations thereof), which, roughly 
speaking, led to a rather weakly or strongly coupled sys¬ 
tems, respectively. 


In the present paper we address this uncertainty 
by utilizing the thermodynamic T -matrix to derive an 
expression for the singlet HQ free energy, Fg^ir). 
Previous work to extract the HQ potential has been 
done both nonperturbatively from IQCD using an effec¬ 
tive Schrodinger equation and potential for the Wilson 
loo p i20l423l] . and within weak-coupling schemes |l8, 
24l25n . Here, we take a diagramatic approach by defin¬ 
ing the potential within a field-theoretic scattering equa¬ 
tion, which leads to a slightly different Schrodinger (or 
Dyson) equation than in previous work. Starting from 
the finite-temperature Bethe-Salpeter equation (BSE) 
for a static quark and antiquark in a mixed coordinate- 
frequency representation, we resum the ladders into a 
closed form expression which naturally serves as the 
definition of the in-medium driving kernel. The per¬ 
tinent spectral function can then be straightforwardly 
related to the HQ free energy. A key element in eval¬ 
uating the resulting expression is the energy depen¬ 
dence of the spectral function, in particular of imagi¬ 
nary parts induced by both one- and two-body correla- 
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tions, i.e., single-quark self-energies and QQ potential- 
type contributions. The former have been calculated be¬ 
fore within the selfconsistent T-matrix formalism ll^ . 
while for the latter we take guidance from perturbative 
calculations 111811241251] . The real part of the potential 
will be taken from a held-theoretic ansatz for a screened 
Cornell potential 1271] . where four parameters charac¬ 
terize its strength and screening properties. This set-up 
will then be deployed to conduct a ht to IQCD data for 
the static color-singlet free energy. 


2. Static QQ Free Energy in Ladder Approximation 


Our objective in this section is to relate the static 
correlation function of (inhnitely massive) quark and 
antiquark to the HQ free energy, FggCr) |i8, 28-33], 
and evaluate it in ladder approximation (BSE or T- 
matrix scheme). Toward this end we hrst formally 
carry out the ladder-diagram resummation of the BSE 
in the imaginary-time (Matsubara) formalism, and then 
reduce the BSE to the T -matrix (potential approxima¬ 
tion). 

The one- part icle propagator in the inhnite-mass limit 
is given by 118 j 


Jn' r' 


,ip-r' 


/ 


d^p' 

(27r)3 iv„ - e, 
1 


1 


(27r)3 h’„ -M- Ed) (v„) 


p,-Ed)(v„,pO 

= b(r')G(')(v„) , 
( 1 ) 


where the b-function implies that the particle never 
moves; E*-') denotes the single-particle selfenergy and 
M is its (heavy) mass; v,j = (2n -H \)nlf5 is a discrete 
fermionic Matsubara frequency (« integer, T - l/y6: 
temperature), adopting the conventions of Ref. iH with 
H - 1. The two-particle Green function inherits the 6- 
function structure M : 

G^(-/r, ri,r 2 |rj,r 2 ) 

= d(ri - r'l) b(r 2 - r') G>(-/t, r ), (2) 

where r - |ri -r 2 |. Here and in the following the “tilde” 
notation is used to denote the reduced Green function 
without the spatial b-functions. The static QQ free en¬ 
ergy, Fqq, can be dehned in terms of the QQ Green 
function as in 

FQQ(r)^-^\n{G>(-i/3,r)) . (3) 


* Color and spin indices are suppressed. 


The task in the following is to calculate G^(-iT, r) in 
Eq. (I2]i within a BSE (or potential model) in the Mat¬ 
subara formalism. In a “mixed” frequency-coordinate 
representation it is given by 


f {Za, V,„ Vm\r) = K (Vn - V,„, 



xG® (za - Vk, Vk\r) f{zA, Vk, 


( 4 ) 


where K (v„ — v^, r) is the interaction kernel and 


GfiZA - Vn, Vn\r) = &^\za - V„)G^^\vn) (5) 


the uncorrelated two-particle propagator. The external 
frequency za is conserved because of imaginary-time 
translation invariance. The full two-particle Green func¬ 
tion takes the form 132. 331 


G{zA,r)^Gf{zA) 

+ ^ ^ Gq\za - Vk, Vk)f(zA, Vk, vi\r)G^Q\zA - v'z, v/), 
P k,i 

(6) 


where we dehned Gf\zA) = -^I,kGf\zA - Vk,Vk), 
and likewise for G(za, Q- The discrete energy depen¬ 
dencies and summations render Eq. (HI of matrix type, 
T„m = K„,n + Y,k[KG^ohnkTkm, which is amenable to 
a formal solution via matrix inversion, T„,n - 2yt[l “ 
KG^Q^]^j^Kk,„, which, however, is not very practical. 

To proceed, we adopt the a potential approximation, 
i.e., neglecting energy transfers in the interaction ker¬ 
nel so that K{v„ - Vk, r) V(za, r)', we keep a pos¬ 
sible dependence on the external energy parameter za, 
in analogy to what has been done, e.g., in the descrip¬ 
tion of electromagnetic plasmas 134] where a “dynami¬ 
cal screening” can induce such a dependence. With this 
approximation, the Matsubara summations in Eq. (HJi 
decouple and one can resum the geometric series to hnd 


f (zA,v„,Vm\r) = 

V(Za, ^ ^k)T{ZA, Vk, Vm\r) 


V(za, r) 


l-V{zA,r)Gf{zA) 


= f{zA\r) 


(7) 


The corresponding full two-particle Green function can 
also be written in compact form as 


G{zA,r) 


1 

[G[f’(zi)]^‘-V(Z2,Q 


( 8 ) 


The above expresson is a key formula of our derivation, 
as it can serve as the dehnition of the interaction kernel 
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in the T -matrix formalism It slightly differs from what 
has been adopted in previous approaches 10 ESS 


25h . as can be seen from inspecting the resulting Dyson 


equation for the meson spectral function. 

Next we use the analyticity of the Green function to 
obtain the correlation function in imaginary-time via a 
spectral representation, 


tion, Eq. (I8]i; we write the pertinent spectral function as 


cry{E, r) - -—Gi{E + 2M + i€\r) » 

7T 

1 -[Vj{E,r) + I.f\E)] 

^[E- Vr(E, r) - + [V,{E, r) + I.f\E)f ’ 

(13) 


G> (-IT, r) = 
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with the spectral function corresponding to Eq. ([8]l, 


a-v(E, r) - —Gj (E -i- ie\r) . (10) 

7T 


where we use to denote the two-particle selfenergy 
that arises from the uncorrelated single-particle selfen¬ 
ergies 1 3 We now make concrete ansatze for the 
potential and selfenergies in terms of a field-theoretic 
model for the in-medium Cornell potential illn an d pre¬ 
vious results from the E-matrix formalism 11261] . The 
former leads to the following real parts of potential and 
uncorrelated selfenergy, 


We use the subscript “V” to indicate that the spectral 
function is calculated in potential approximation. Here 
and in the following, the subscripts ”R” and are 
short-hand notations for the real and imaginary part of 
a given quantity. Since we work in the static limit 
(M —> oo), we redefine the spectral function relative 
to the QQ threshold by a shift E - E' - 2M, i.e., 
&v{E, r) = cry(E + 2M, r), which further simplifies 
Eq. @ to 


G^(-iT, r) 


CO 

9-L^ ^ ^ dEcTy(E,r)Q-^^ ( 11 ) 

— CO 


Combining Eqs. Q, (fTOl i and (fTTT i. we arrive at the final 
expression for FqqG) within the T -matrix approach. 


PQQ = -p In 


CO 

/ 


dE&y(E, 


( 12 ) 


In the following section, in the spirit of a variational 
approach, we will fit a physically motivated ansatz for 
the potential to IQCD data of the in-medium free energy. 


3. Potential Extraction from the Free Energy 

The two-particle Green function in medium receives 
contributions from the (medium-modified) interaction 
kernel, V, and from the single-particle selfenergies, 
arising from individual interactions of each particle with 
the heat bath which figure into Gf\ Both are, in prin¬ 
ciple, complex quantities and are contained in our T- 
matrix expression for the full two-particle Green func- 


4 

VR(E,r) = --a, - 

3 r 


Z®(£) = 


■-ocsmo + a -— 
3 m, 


(14) 


(15) 


corresponding to the r-dependent and r-independent 
terms of the screened Cornell potential, respectively 
(additional contributions from the uncorrelated 2- 
particle propagator are small and therefore ne¬ 
glected in the present exploratory calculation). Their 
strengths are characterized by the strong coupling, as - 
I (An), and string tension cr, along with correspond¬ 
ing screening masses nio and m,, respectively. As in 
our previous work il2ll . we utilize the latter as indepen¬ 
dent fit parameters, where, for simplicity, we choose a 
linear ansatze, nio s - CD,sgT with two constant coef¬ 
ficients, cd.s (which turn out to be close of order one). 
The imaginary part of the uncorrelated 2-particle self- 
energy, Z} , follows from the folding integral of the 
single-particle propagators which contain the respec¬ 
tive in-medium selfenergies. In principle, this quan¬ 
tity follows automatically from a selfconsistent calcu¬ 
lation, once the potential is specified, as was carried out 
in Ref. ll26ll . For our exploratory study in the present 
work, we do not enforce selfconsistency but instead em¬ 
ploy the energy dependence of Z^*^ as found in Ref. 

(see Fig. 3 in there) to compute \ It turns out that the 
latter can be well represented by Gaussian form. 


Z®(E) = E>exp 


I® 

(E - ^)2 


2(CE)2 


4 

-OsT . 


(16) 


With C ^ 1.6 we can reproduce the energy dependence 
of the single-particle width of Ref. being very sim¬ 
ilar for both internal and free energies as underlying po¬ 
tential. The strength parameter D, however, which we 
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define relative to the perturbative value, ^OsT, markedly 
depends on the underlying potential, varying by a factor 
of 2-3 between free and internal energies; as mentioned 
above, we here use it as a fit parameter. 

Finally, we have to specify the imaginary part of the 
potential; its magnitude and r dependence has been 
determined previously in the weak coupling limit as 


= %asTQ>(mDr) ifl^ 1^ where <[>(x) is a 


smoothly varying function which vanishes for x=0 and 
approaches 1 for large x. However, in our definition of 
the potential, we have identified the non-zero part at in¬ 
finite distance with the uncorrelated selfenergy. Upon 
subtracting this part and again allowing for a nonpertur- 
bative enhancement with the same stength parameter D 
and energy dependence as in Eq. (fTbl l we arrive at the 
form 


Vi{E, r) - Dexp 


{E - 
2{CTf 


-aJ{Q>{BmDr) - 1 ] 


(17) 

To accommodate the possibility that the imaginary part 
of the string interaction can have a different r depen¬ 
dence than given by <1) (albeit with the same asymptotic 
limits), we rescaled its argument by a constant B (which 
turns out to be close to one). We note that nontrivial en- 
ergy dependencies in both imaginary parts, V; and , 
are crucial to ensure convergence of the integration in 
Ea. dTSt . In this sense, the potential is required to go 
beyond the static approximation. 

As an example of our method we perform a fit to 
a set of free energies computed in lattice QCD with 
Nf - 2 + I dynamical flavors 1^ . for temperatures 
T^l.2, 1.7, 2.2 Tc, with 7’c=196MeV Q. The fol¬ 
lowing parameter values were used; The Coulomb De¬ 
bye mass is nio - 1 ■22gT, close to the expected HTL 
value, while the string interaction is screened less with 
nis - O.SgT; for the string tension itself a constant value 
of cr = 1.58GeV/fm turns out to be sufficient, which 
is larger than in vacuum but significantly smaller than 
in previous fits to the in-medium free energy lll2l UTT . 
Eor the strong coupling we utilize a weak temperature 
running, = 0.288 - 0.05 (yt^ - 1), decreasing from 
0.288 at 1.27’c to 0.238 at 2.4 T^. Eor the dimensionless 
coefficients we have B = 0.73, which delays the ap¬ 
proach of O to its asymptotic value of one, and C-2.Q7 
in Eq. (fTbl l. which produces slightly longer tails in the 
energy dependence of the HQ selfenergies than in the 
microscopic calculations of Ref. lE^ . 

The resulting fit to the IQCD free energies is dis¬ 
played in Eig.[Tl showing good agreement. Reproducing 
the temperature dependence essentially ensures that the 


internal energies, U = F + TS with S - —dFjdT, also 
agree with the lattice data. The underlying two-body 
potential turns out to lie significantly above Fqq, but 
below Uqq, at intermediate and long distances. How¬ 
ever, the potential differs from both free and internal en¬ 
ergy in that its derivative, characterizing the inter-quark 
force, is larger than for both Fqq and Uqq, in a region 
around r ^ 1 fm. This is a remnant of the confining 
string term in the hot medium, caused by the relatively 
small screening mass, nis- Even though this force is not 
so large, it has a rather long range, encompassing a rela¬ 
tively large volume. In the heavy-light quark sector, this 
enables heavy quarks to effectively couple to more parti¬ 
cles in the heat bath, relative to a short-range force. Our 
preliminary calculations of the HQ diffusion coefficient 
indeed indicate that the latter is comparable to earlier 
calculations using the internal energy ifl^ . which has a 
stronger force but over a shorter range. In this sense, the 
potential obtained in the present study may still imply a 
strongly coupled QGR 

The above exmaple of a fit is not unique, mostly due 
to uncertainties in the imaginary parts, which in the fu¬ 
ture will have to be calculated in a selfconsistent scheme 
as laid out in Ref. ME However, several qualitative, 
model-independent observations can be made. Eirst, if 
the imaginary parts, iV + iP'^)j(F), are taken to zero, 
the spectral function cry becomes a 6 function in energy 
and the potential in the denominator is forced to become 
the free energy, thus recovering the weak-coupling anal¬ 
ysis of Ref. il8ll . This underlines the key role of the 
imaginary parts of both potential and single-quark self¬ 
energy in raising the potential above the free energy, 
i.e., rendering a stronger force. The E-depedence of 
(V + iP'^)i(E) (the “width of the width”), encoded in the 
width parameter C, also plays an important role, which 
reiterates the need for microscopic calculations, such as 
the T -matrix approach. The dependence on the real part 
is more robust and “model independent” in the sense 
that it can be largely determined by the small distance 
behavior of free energy where <h(x) forces the imagi¬ 
nary part to become small. Let us also comment on the 
magnitude of the string tension exceeding its vacuum 
value by about 60%. This appears to be a rather generic 
feature within the screened Cornell potential ansatz. In 
previous works utilizing this ansatz lll2ll27n . where the 
strength of the string term was additionally allowed to 
develop a temperature dependence, fits to the free en¬ 
ergy preferred string tensions of up to a factor 3-4 times 


^The imaginary parts implicit in our fit are comparable to what was 
found before with the internal energy as potential, which indicates that 
they are in the expected ballpark of the selfconsistent solution. 
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Figure 1: Our fit (blue lines) to the color-singlet QQ free energies computed in lattice QCD (black dots) 1^ at 3 different temperatures, 7= 1.2, 1.7, 
2.2 Tc in the left, middle and right panel, respectively. The red lines are the real parts of the underlying potential, while the green lines represent 
the correpsonding internal energies. 


the vacuum, albeit with a much stronger screening than 
in our present potential fit. As a result, the in-medium 
free energy never exceeds the vacuum potential (as in 
our present work). Similar features are also found in ap¬ 
plications of the Dyson-Schwinger approach to quarko- 
nia at finite T 1^ . A possible microscopic mechanism 
for an enhanced string tension could be the release of 
magnetic monopole in the near-Tc region, as suggested, 
e.g., in Ref. 1^ . 

4. Conclusions and Outlook 

Utilizing a thermodynamic T-matrix approach, we 
have derived an expression for the static heavy-quark 
free energy that directly relates to the underlying in¬ 
teraction kernel, resummed in ladder approximation. 
The such defined kernel serves as our definition of an 
in-medium potential, VqqCt), in the context of a dia¬ 
grammatic many-body approach. Key elements of this 
framework are a non-trivial energy dependence of the 
pertinent QQ spectral function, induced by the imagi¬ 
nary parts of both one- and two-body type, i.e., in the 
single-quark selfenergies and the potential. We have 
applied this set-up to fit lattice-QCD “data” for the 
color-singlet free energy in the spirit of a variational 
scheme for the underlying potential. To facilitate this, 
we made an ansatz for the real part of the potential in 
terms of a screened Cornell potential, augmented with 
an imaginary part motivated by perturbative studies but 
with variable magnitude. We additionally accounted for 
complex heavy-quark self-energies with an energy de¬ 
pendence motivated by earlier selfconsistent implemen¬ 
tations of the many-body theory. The resummations 
of the interaction kernel, with a nonperturbative string 
term, and of the selfenergies are required in a regime of 
strong coupling, as expected for the QGP at tempera¬ 
tures not too far above Tc- Indeed, our resulting poten¬ 
tial significantly deviates from the weak coupling limit 


(in which case it would be close to the free energy), and 
exhibits a rather long-range remnant of the confining 
force surviving well above Tc- The associated imagi¬ 
nary parts are large, not inconsistent with what was ob¬ 
tained earlier using the internal energy as a potential ap¬ 
proximation. First estimates indicate that the pertinent 
HQ diffusion coefficient is around Ds{2nT) ^ 3 - 5 , 
suggestive for a strong coupling regime. 

Several further investigations are in order to scruti¬ 
nize our initial estimates within the T -matrix approach. 
First, the calculation needs to be carried out selfconsis- 
tently, by computing the heavy-quark selfenergies from 
the underlying T -matrices and reinserting the former in 
the latter 11261] . Second, the free energies should be up¬ 
dated with the most recent lattice-QCD results. Third, 
quarkonium correlators, HQ suscpetibilities and trans¬ 
port coefficients should be computed and systematically 
compared to lattice data. The role of possible three- 
body correlations may need to be addressed. If the the¬ 
oretical framework passes these tests, applications to 
heavy-quark and quarkonium phenomenology are war¬ 
ranted. 
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